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X 



D.J. Singh and I.I. Mazin 
Center for Computational Materials Science, 
Naval Research Laboratory, Washington, DC 20375, U.S.A. 
(February 1, 2008) 

The electronic structure, stability, electron phonon coupling and superconductivity of the non- 
oxide perovskite MgCNi3 are studied using density functional calculations. The band structure is 
dominated by a Ni d derived density of states peak just below the Fermi energy, which leads to a 
moderate Stoner enhancement, placing MgCNi3 in the range where spin fluctuations may noticeably 
affect transport, specific heat and superconductivity, providing a mechanism for reconciling various 
measures of the coupling A. Strong electron phonon interactions are found for the octahedral rotation 
mode and may exist for other bond angle bending modes. The Fermi surface contains nearly 
cancelling hole and electron sheets that give unusual behavior of transport quantities particularly 
the thermopower. The results are discussed in relation to the superconductivity of MgCNi3. 



The discovery^ of superconductivity in the non-oxide 
perovskite MgCNi3, with critical temperature T c «8K, 
raises the questions of how superconductivity appears in 
such a Ni rich phase and its relationship to other unusual 
superconducting phases, partkjularly the NL, Pd and Pt 
borocarbides and boronitridestm and PdHD u These pos- 
sibly related materials have relatively high values of T c 
and high concentrations of magnetic or near magnetic 
group 8 elements but different underlying physics. 

Typically, but not always, perovskites are distorted 
by the freezing in of unstable zone boundary rotational 
phonons to form structures like the Puma, GdFeC>3 
structure or off-centerings like those that produce ferro- 
electricity as in BaTi03. This is due to the importance 
of ionic and central interactions in many perovskites; the 
stability of the ideal cubic structure depends on a balance 
between A-site cation-anion and the i?-site cation-anion 
interactions. In the minority that are truely cubic, there 
are often soft anharmonic phonons of displacive charac- 
ter, as e.g. in KTaC>3 and/or rotational character (the 
combination common in Pb based perovskites). 

According to refinements ,0 MgCNi3 is stoichiometric 
with only a very small, 4% carbon deficiency and occurs 
in the ideal cubic perovskite structure. The ambient tem- 
perature lattice parameter is 3.812A, which yields a C-Ni 
bond length of 1.906A. Considering that Mg is expected 
to occur as Mg 2+ and so the CM3 subunit as negatively 
charged, this bond may be considered short. As such, 
strong Ni-C covalent interactions are expected. If so 
and if bonding or antibonding electronic states associ- 
ated with this are present at the Fermi energy Ep, strong 
electron phonon coupling (EPC) may be expected. The- 
oretical studies of the electronic structure and phonons 
showed that in fact EPCs associated with relatively hard 
modes that modulate strong bonds involving Ni and first 
row elements play an important role inpttie-superconduc- 
tivity of the borocarbides and nitrides.EJHJ 

On the other hand, the perovskite topology and the 
fact that Mg likely acts as a simple cation that is elec- 
tronically inactive near Ep suggests that perhaps low 
frequency anharmonic modes play an important role. In 



particular, phonon modes, like that associated with the 
octahedral rotation, will modulate the C-Ni-C bond an- 
gles, which in this structure are expected to be impor- 
tant in the hopping associated with metallic conduction. 
The crystal structure may be conceptually viewed as ex- 
panded fee Ni with 25% of the sites replaced by Mg 
and C interstitials in octahedral sites. The bands thus 
should be narrow and transition metal-like around Ep 
with a higher filling than in pure Ni. This, suggests some 
analogy with the band structure of PdHM, also an ex- 
panded group 8 metal with a higher filling. The main 
difference is the importnat role of C as hopping mediator 
in MgCNi 3 , as discussed below. 

The present experimental situation is confusing. Tun- 
neling and upper critical field measurements were re- 
ported by Mao et al£3 These characterize MgCNi3 as a 
strong coupling superconductor with an unusually large 
reduced enerjpz.gap, 2A(0)//cT c « 10, near the upper the- 
oretical limitlij for s-wave superconductivity. However, 
the specific heat junan AC(T c )/kT c w 1.9 suggests mod- 
erate A <1 couplingJlil Furthermore, they suggested that 
MgNiC3 may be a non-s superconductor, based on ob- 
servation of a zero-bias anomaly (often, but not always, 
due to a sign-changing order parameter) . This implies an 
analogy with metallic Pd where superconductivity is not 
observed, but a possible spin fluctuation mediated mech- 
anism was discussed, or Sr2Ru04, which is apparently a 
triplet superconductor due to spin fluctuations .113 On the 
other hand, doping with Co, a likely magnetic impurity, 
rapidly decreases the superconducting volume Jij'O, and 
the same effect, but stronger, is found with MnO Cu 
doping, however, reduces, T c without suppressing the su- 
perconducting fraction.E3 This is consistent with expec- 
tations for -s.- wave pairing in a rigid band model. Hay- 
ward et alS3 reported a calculated electronic density of 
states (DOS) that shows a large peak just below Ep 
and speculated about tie importance of magnetic fluc- 
tuations based on this.E3 In this regard, doping on the 
Mg j4-site may be illuminating as rigid band behavior 
may be more likely in that case. 

Here we report electronic structure calculations focus- 
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ing on these issues. These were done in the local density 
approximation (LDA)t3 with the general potential lin- 
earized augmented planewave (LAPW) method.Ej Well 
converged gone samplings and basis sets, including lo- 
cal orbitalsE] to treat the Mg semicore states and relax 
linearization errors for the Ni d bands were used. The 
experimental lattice constant of a=3.82A was used (our 
calculated LDA lattice parameter is 1.7% smaller as is 
typical for this approximation) . We study the band struc- 
ture and Fermiology, the susceptibility and proximity to 
magnetism, transport coefficients and selected phonon 
modes and electron-phonon couplings. Transport coef- 
ficients were determined using kinetic theorjO with zone 
integrations using LAPW eigenvalues on a grid of 2925 
k-points in the irreducible 1/48 wedge. 

The band structure and DOS are shown in Figs. ^ and 
0. The Fermi surfaces are in Fig. ||. As expected, Mg 
plays an minor role in the bands in the valence region. 
(There is a weak Ni d-Mg s bonding interaction. The T 
point state at -6 eV has Mg-Ni bonding character with 
the corresponding antibonding state at +3 eV. This co- 
valency is localized around T). Octahedral coordination 
is unfavorable for formation of C sp hybrids, so not sur- 
prisingly the C s orbitals are also inactive. They occur 
around -12 eV (relative to Ep) and are not shown. The 
valence region, which extends from -7.1 eV to +1 eV is 
therefore derived from the 15 Ni d and 3 C p orbitals filled 
by 34 electrons. The C p orbitals are strongly hybridized 
with Ni d, and are located below most of the d states. 
So in the first approximation they can be integrated out 
and the bands near the Fermi level can be analyzed in 
terms of the Ni d states. 

The unusual two-fold linear Ni coordination by C 
makes some of the bands very narrow. For instance, 
Ni(x) yz and y 2 — z 2 , Ni(y) zx and z 2 — x 2 , and Ni(x) xy 
and x 2 — y 2 orbitals do not disperse in the nearest neigh- 
bor approximation (NNA) as they have no C orbitals to 
hop through.Ea The remaining 9 orbitals form three inde- 
pendent manifolds, consisting of Ni(x) xy, Ni(z) yz, and 
Ni(y) y 2 — r 2 orbitals, coupled via C p y and the two corre- 
sponding combinations (i.e. p x andp 2 ). Interestingly, in 
the NNA, one of the three resulting bands in each man- 
ifold, the one that involves the antibonding combination 
of the t2 9 orbitals, is non-bonding again. The result is 3 
bonding, 3 nonbonding and 3 antibonding bands. Two 
of these antibonding bands cross Ep. In the large crystal 
field limit (relative to the small widths), these have a sim- 
ple dispersion proportional to svn(ak x /2) 2 + sm(ak y /2) 2 
etc., which naturally makes them flat for a number of 
high-symmetry directions. This is reflected in the Fermi 
surface (Fig. ||). The lower band produces T centered 
rounded cube shaped electron sections and a narrow low 
weight jungle gym along the R — M lines also contain- 
ing electrons. The upper band (Fig. || bottom) forms 
a two-sheet Fermi surface consisting of pancaked squares 
centered at the X points and ovoids along the T — R lines, 
both hole-like. The flat square shape reflects weak dis- 
persion of the underlying band along X — M and strong 



dispersion along X — T (The quasi-2D behavior is be- 
cause sm(ak x /2) 2 + sm(ak y /2) 2 does not disperse along 
z.) Near T, where these bands make up the three fold 
degenerate state (this is the second one below Ep), the 
dispersion is due to direct Ni-Ni dda hopping. Turning 
to the two bands that cross Ep, the lower band forms 
the more interesting heavy sheet of Fermi surface around 
X. At an M point, say (110), sm(ak x /2) 2 + sm(ak z /2) 2 
and sin(ak y /2) 2 + sin(afc 2 /2) 2 are degenerate and their 
hybridization via Ni(x) xy and Ni(y) xy orbitals is pro- 
portional to cos(ak x /2) cos(ak y /2) and vanishes near M, 
which is why they are heavy. On lowering Ep, the X cen- 
tered squares grow grow, as does DOS at Ep, N(Ep), un- 
til they meet at the M points and the topology changes. 

The value of N(Ef) is 4.99 eV -1 on a per formula 
unit basis. Two recent non- full-potential studies ob- 
tain differing values. Dugdale and Jarlborg report 6.35 
eV -1 and 3.49 eV -1 depending on the exact method 
they employcJ while Shim and Min obtaift,5.34 eV -1 , 
which is much closer to the present result .Ej These dif- 
ferences are significant as they control the proximity to 
magnetism. In particular, the susceptibility, x(0), is de- 
termined by N(Ep)/(l — N(Ep)I), where / is the Stoner 
parameter, and the denominator contains a small differ- 
ence involving N(Ep). Comparing our calculated full 
potential N(Ep) with the experimental linear coefficient 
7 «10 mJ/moleNiK 2 , we obtain a specific heat renormal- 
ization 7 /j band =2.6. The calculated plasma frequency is 
hujp = 3.25 eV. The calculated Hall number is -1.3 xlO 22 
cnu. 3 . This agrees well with the measured value of Li et 
al.cB The constant scattering time approximation ther- 
mopower, S, is p-type except at very low T below 10K 
where it is n-type consistent with the Hall number. S is 
very small (less than 1 fiX/K) below 150K but then rises 
more rapidly reaching 5 fiV/K at 300K and 16 pV/K at 
600K. Both the Hall number and the thermopower are 
controlled by the competition between the hole and elec- 
tron pockets of comparable size, which leads, for instance, 
to the unusual T— dependence of thermopowpc— Using 
the calculated top and the measured resistivityaO in the 
Bloch-Gruneisen formula, dp/dT=(8ir 2 /hujp)kB^tr, we 
obtain X tr ~ 1 — 1.6, though we note that this extraction 
of Xtr is sensitive to sample quality. 

The simplest estimate of A from zone center frozen- 
phonon calculations is not applicable to cubic perovskites 
since no zone-center modes couple by symmetry. How- 
ever one can use frozen-phonon zone-corner calculations 
for this purpose, as in Ref. ^(| The relevant formula is 

1 V 1" Lu 'j/q 

1 v n,m, k,q 

= 27V T £(|.gf M, q > 

Here is the per spin DOS at Ep, and g v is an 
electron-ion matrix element from the derivative of the 
ionic potential with respect to the dimensionless phonon 
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coordinate (see Ref. g6|for details). To get qualitative in- 
formation on the EPC, we focus on two R-point phonons. 
These are the octahedral rotation and the fully symmet- 
ric breathing mode. The former mode changes the C- 
Ni-C bond angles, but not the bond lengths in lowest 
order, while the latter is a pure bond stretching mode. 
The calculated frequencies are 105 cm -1 and 349 cm -1 , 
respectively. 

To estimate the matrix element, we selected several 
points on the intersection line of the k and k+q Fermi 
surfaces, and fitted the bands in the nearby region with 
the second-order perturbation theory, and then averaged 
the resulting matrix element. We included three inter- 
section points, where g does not vanish by symmetry, 
for the breathing mode, and 7 points for the rotational 
mode. Assuming that averaging over q does not change 
g 2 , we find for the breathing mode A » 0.005 and for the 
rotational mode A ss 1.2. The latter dominates not only 
because of its nearly 4 times larger deformation potential 
and its 15 times smaller denominator Mio 2 . Because of 
the mode degeneracy, this corresponds to a total rota- 
tional modes X ro t ~ 3.6 However, in oxide perovskites, 
this i?25 rotation mode usually stiffens rapidly away from 
the zone boundary, reflecting the rigidity of the O octahe- 
dra. Here such a stiffening may also be expected. Thus, 
the zone averaged A will likely be considerable smaller, 
but probably still substantial. In any case, it can be said 
that the stiff C-Ni bond stretching modes are apparently 
not significant contributors to the EPC, while rotational 
(and probably other C-Ni-C bond bending) modes are 
strongly coupled. Note that the popular rigid muffin-tin 
potential method is hard to apply here: for transition 
metals, it is known to overestimate EPC, and Ni is in a 
low-symmetry position. Most importantly, with such a 
large disparity of the contributions from different modes, 
there is no telling beforehand, which average phonon fre- 
quency should be used in calculating A. 

As mentioned, possible nearness a magnetism could 
conceivably play an important role. To investigate this, 
we have calculated spin susceptibility directly from the 
variation of the total energy with small imposed mag- 
netizations (ranging from 0.1 to 0.7 (1b per cell) in 
fixed spin moment calculations, SE to t = x^m 2 /2 = 
(TV^T 1 — I)m 2 /A, where / is the Stoner factor, characteriz- 
ing intraatomic exchange (for a compound like MgCNi3, 
I is expected to be close to 1/3 of the pure Ni value of 
1.16 eV, probably slightly reduced because of hybridiza- 
tion with C). It appears that 8E to t{m) noticeably de- 
viates from the quadratic behavior. This is expected in 
extended Stoner theoryJ23 where an average N(m) is sub- 
stituted for Nf. With N{m), the above expression gives 
a good description of SE tot (m) up to m < 0.6 fiB with a 
weakly energy dependent / w 0.95/3 eV. Self-consistent 
virtual crystal calculations corresponding to a 10% re- 
placement of Mg by a monovalent element (Na in the 
calculations but with the lattice parameter fixed - so Li 
or a 5% Mg deficiency may be a better experimental 



analogue) produced a borderline ferromagnetic ground 
state with a very low energy gain and a moment of 0.25 
/is/cell. 

The susceptibility renormalization (1 — IN^)^ 1 » 5, 
and \ — 2.7 x 10~ 4 emu/moleNi_ somewhat larger 
than the reported experimental valueEj of 1.7/timeslO -4 
emu/moleNi. This number is expected to be very sensi- 
tive to doping, i.e. composition, due to the near cance- 
lation in the denominator. Perhaps the difference with 
experiment can be understood in these terms e.g. if the 
4% of C vacancies lead to a higher effective band fill- 
ing. If the Stoner renormalization is indeed rj 5, this 
signals presence of significant spin fluctuations. To com- 
pare with, in Sr2Ru04, where (1 — IN^)~ X rj 9, spin fluc- 
tuations are believed to cause triplet superconductivity, 
and in Pd metal, where it is « 30, spin fluctuations de- 
stroy superconductivity which would otherwise exist due 
to the sizeable (A > 0.5) EPC. Importantly, the effect of 
spin fluctuations on mass renormalization, transport and 
superconducting properties is very different: essentially, 
in the first two cases, the coupling constants add, while 
for superconductivity they add in the mass renormaliza- 
tion term, and subtract in the pairing term (In the strong 
coupling T c is not proportional to exp[— (1 + A) /A], but 

rather GXp [ (l ~h ^ph H~ ^spin ) / \^-ph ^spin 

)].) This may 

explain the inconsistency between coupling constants de- 
termined from different experiments. 

In any case, the present results underscore certain sim- 
ilarities with PdH in terms of band structure and mag- 
netic renormalizations, as well as with the borocarbides. 
However, the fact that the most important phonons are 
Ni bond-bending modes places MgCNi3 in a unique class 
of its own. 
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FIG. 2. Electronic DOS (upper panel) of MgCNi 3 as 
calculated within the LDA. The dashed line is the d con- 
tribution within the Ni LAPW spheres of radius 2.04 ao. 
The lower panel is a blow-up around Ef- 
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FIG. 1. LDA band structure of MgCNi 3 . E F is at 0. 
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FIG. 3. Calculated Fermi surfaces. The top panel is 
from the lower band and shows eight zones. The rounded 
cube sections are centered at F, while the thin jungle gym 
is along the M — R lines. The bottom panel is from the 
upper band, showing dimpled square shaped sections cen- 
tered at X and ovoids along F — R. 
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